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Abstract 

A universal period doubling cascade analogous to the famous Feigenbaum- 
Coullet-Tresser period doubling has been observed in area-preserving maps 
of M. 2 . Existence of the "universal" map with orbits of all binary periods 
has been proved via a renormalization approach in (Eckmann et al 1984) 
and (Gaidashev et al 2011). These proofs use "hard" computer assistance. 

In this paper we attempt to reduce computer assistance in the argument, 
and present a mild computer aided proof of the analyticity and compact- 
ness of the renormalization operator in a neighborhood of a renormalization 
fixed point: that is a proof that does not use generalizations of interval 
arithmetics to functional spaces — but rather relies on interval arithmetics 
on real numbers only to estimate otherwise explicit expressions. The proof 
relies on several instance of the Contraction Mapping Principle, which is, 
again, verified via mild computer assistance. 
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1 Introduction: Period doubling universality for area- 
preserving maps 

Following the pioneering discovery of the Feigenbaum-Coullet-Tresser period dou- 
bling universality in unimodal maps (Feigenbaum 1978, Feigenbaum 1979, Tresser 
and Coullet 1978) universality has been demonstrated to be a rather generic phe- 
nomenon in dynamics. 

To prove universality one usually introduces a renormalization operator in 
a functional space, and demonstrates that this operator has a hyperbolic fixed 
point. 

Such renormalization approach to universality has been very successful in 
one-dimensional dynamics, and has led to explanation of universality in unimodal 
maps (Epstein 1986, Epstein 1989, Sullivan 1992, Lyubich 1999), critical circle 
maps (de Faria 1992, de Faria 1999, Yampolsky 2002, Yampolsky 2003) and holo- 
morphic maps with a Siegel disks (McMullen 1998, Yampolsky 2007, Gaidashev 
and Yampolsky 2007) . Universality has been also abundantly observed in higher 
dimensions, in particular, in two and more dimensional dissipative systems (cf. 
(Collet et al 1980), (Sparrow 1982)), in area-preserving maps, both as the pe- 
riod doubling universality (Derrida and Pomeau 1980, Helleman 1980, Benettin et 
al 1980, Collet et al 1981, Eckmann et al 1982, Eckmann et al 1984, Gaidashev and 
Koch 2008, Gaidashev and Johnson 2009, Gaidashev and Johnson 2009, Gaida- 
shev et al 2011), and as the universality associated with the break-up of invariant 
surfaces (Shenker and Kadanoff 1982, MacKay 1982, MacKay 1983, Mehr and 
Escande 1984), and in Hamiltonian flows (Escande and Doveil 1981, Abad et 
al 2000, Abad et al 1998, Koch 2002, Koch 2004, Koch 2008, Gaidashev and 
Koch 2004, Gaidashev 2005, Kocic 2005, Khanin et al 2007). 

In this paper we will consider a period doubling universality for area-preserving 
maps of the plane — an analogue of Feigenbaum-Coullet-Tresser universality in 
higher dimensions. 

An infinite period doubling cascade in families of area-preserving maps was 
observed by several authors in early 80's (Derrida and Pomeau 1980, Helleman 
1980, Benettin et al 1980, Bountis 1981, Collet et al 1981). A typical period 
doubling scenario can be illustrated with the area-preserving Henon family (cf. 
(Bountis 1981)) : 

H a (x,y) = {-y + 1 - ax 2 ,x). 



Maps in this family posses a fixed point ((—1 + a/1 + o)/a, (—1 + vT+~a)/a) 
which is stable for —1 < a < 3. When a\ = 3 this fixed point becomes unstable, 
at the same time an orbit of period two is born with H a (x±, x^) = (x^,x±), 
x± = (1 ± \/a — 3) /a. This orbit, in turn, becomes unstable at 02 = 4, giving 
birth to a period 4 stable orbit. Generally, there exists a sequence of parameter 
values a*;, at which the orbit of period 2 k ~ 1 turns unstable, while at the same 
time a stable orbit of period 2 k is born. The parameter values accumulate on 
some Oqo. The crucial observation is that the accumulation rate 
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is universal for a large class of families, not necessarily Henon. 

Furthermore, the 2 k periodic orbits scale asymptotically with two scaling 
parameters 

A = -0.249..., n = 0.061. .. 

To explain how orbits scale with A and \i we will follow (Bountis 1981). 
Consider an interval (afc,afc+i) of parameter values in a '"typical" family F a . For 
any value a € (dk,o-k+i) the map F a posses a stable periodic orbit of period 
2 k . We fix some otk within the interval (afc,cifc+i) in some consistent way; for 
instance, by requiring that the restriction of F* k to a neighborhood of a stable 
periodic point in the 2 fe -periodic orbit is conjugate via a diffeomorphism H k to 
a rotation with some fixed rotation number r. Let p' k be some unstable periodic 
point in the 2 fc_1 -periodic orbit, and let p k be the further (from p' k ) of the two 
stable 2 fc -periodic points that bifurcated from p' k . Then, 

I = - Mm |pfc ~^' * = - lim 

A k^oo \pk+l - P' k+1 \ jU k^oo p k+ i 

where p k is the ratio of the eigenvalues of DH k {p k ). 

This universality can be explained rigorously if one shows that the renormal- 
ization operator 

R[F] = A" 1 0F0F0 A F , (1.1) 

where is some F-dependent coordinate transformation, has a fixed point, and 
the derivative of this operator is hyperbolic at this fixed point. 

It has been argued in (Collet et al 1981) that A^ is a diagonal linear trans- 
formation. Furthermore, such Ap has been used in (Eckmann et al 1982) and 
(Eckmann et al 1984) in a computer assisted proof of existence of the renormal- 
ization fixed point. This is the strongest result to date concerning the existence 
of the renormalization fixed point. 

The problem has been further studied in (Gaidashev and Koch 2008), where 
the authors have demonstrated that the period doubling renormalization fixed 
point of (Eckmann et al 1984) is "almost" one dimensional, in the sense that it 
is close to the following Henon-like map: 

H*(x, u) = {4>(x) — u, x — (f>(4>{x) — u)), 

where <\> solves 

2 

4>( x ) = jHH^x)) ~ x - 

(Gaidashev and Koch 2008) also suggests an analytic approach to the full 
period doubling problem for area-preserving maps based on its proximity to the 
one dimensional. 

Infinitely renormalizable maps, that is maps on the renormalization stable 
manifold, have been studied in (Gaidashev and Johnson 2009) and (Gaidashev 
et al 2011), where it has been shown that such maps (in a codimension 1 stable 
submanifold) posses stable Cantor sets (sets on which the Lyapunov exponents 
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are equal to zero) on which dynamics is rigid: the stable dynamics for two differ- 
ent infinitely renormalizable maps is conjugated by a C 1+a transformation (that 
is, a differentiable transformation whose derivative is Holder with exponent a). 
This is in a stark contrast to the situation in dissipative Henon-like maps, where 
universality of Cantor sets (around the so called "tip" ) coexists with non-rigidity 
of the conjugacy of the dynamics on these sets (see (de Carvalho et al 2005, Lyu- 
bich and Martens 2008, Hazard et al 2010)). 

The computer-assisted proof of (Eckmann et al 1984) can be described as 
an application of the Contraction Mapping Principle for a Newton map for the 
operator (|l.ip in some neighborhood in an appropriately chosen Banach space. 
The proof is similar in spirit to the proofs (Lanford 1982, Lanford 1984) of the 
Feigenbaum-Coullet-Tresser universality by Oscar Lanford, and uses bounds for 
analytic functions (not necessarily polynomials) via the so called "standard sets" : 
loosely, generalizations of intervals in K to infinite-dimensional functional spaces. 
Such computer assisted proofs came to be called "hard" within the computa- 
tional community. The approach turned out to be quite powerful, albeit rather 
specialized and time and effort consuming. 

In this paper we study analyticity and compactness of the renormalization 
operator (jl.ip using, almost exclusively, the interval arithmetics in M. We ob- 
tain a bound on the neighborhood of analyticity of the operator. Specifically, 
we identify a low order polynomial approximation for the "generating function" 
of the fixed point map of (jl.ip . and estimate the radius of the neighborhood of 
analyticity of the renormalization operator, centered on this polynomial approxi- 
mation in an appropriate Banach space of functions analytic on a certain domain 



Furthermore, we demonstrate that renormalizations of functions from this 
neighborhood of analyticity, are themselves functions, analytic on a larger domain 
in M 2 , which implies compactness of the renormalization operator. 

We achieve these two goals using only what can be termed as "mild" computer 
assistance: computer estimates on arithmetic operations on reals and on the 
square root function on reals, as well as norm estimates of low order (quadratic) 
polynomials. It is also our hope that the reader can find of interest several little 
"tricks" used in the paper to get bounds on solutions of functional equations. 

2 A renormalization operator on generating functions 

An "area-preserving map" will mean an exact symplectic diffeomorphism of a 
subset of M 2 onto its image. 

Recall, that an area-preserving map F : (x, u) t-t (y, v) can be uniquely 
specified by its generating function S: 




(2.2) 



if the equation u 



Si(x,y) has the unique solution y(x,u). 



1 



Furthermore, we will assume that F is reversible, that is 

ToFoT = F~ 1 , where T(x, u) = (x, — u). 
For such maps it follows from (|2.2[) that 

Si(y,x) = S 2 {x,y) = s(x,y), 

and 

/ x \ F ( y 
\-s(y,x)J h-> \s(x,y) j 

It is this "little" s that will be referred to below as "the generating function" . 
Applying a reversible F twice we get 

x' \ F f z> \ = f z> \ F f y> \ 
-s(z',x')J ' y \s(x',z')J \-s(y',z')J i y \s(z',y')J ' 

It has been argued in (Collet et al 1981) that 

Ap(x,u) = (Xx,fiu), 

where A and \x are some real parameters, dependent on F. 

We therefore set (x',y') = (Ax, Ay), z'(Xx,Xy) = z(x,y) to obtain: 



( x \ A ( Xx \ FoF f x y \ A ^( y \ r2 on 

\-±s(z,\x)J h> \-s(z,\x)J ' y \s(z,Xy)J ' y [j l s(z,Xy)J , K '' 
where z(x, y) solves 

s(Xx, z(x, y)) + s(Xy, z(x, y)) = 0. (2.4) 

If the solution of (j2.4[) is unique, then z(x,y) = z(y,x), and it follows from 
(|2.3p that the generating function of the renormalized F is given by 

s(x, y) = n~ 1 s(z(x, y), Xy). (2.5) 

Furthermore, it is possible to fix some normalization conditions for s and z 
which serve to determine scalings A and \x as functions of s. Notice, that the 
normalization 

s(l,0) = 

is reproduced for s as long as 

*(l,0)=z(0,l) = l. 
In particular, this implies that 

s(A,l) + s(0,l) = 0. 



2 



Furthermore, the condition 

Mi>o) = i 

is reproduced as long as 

H = d lZ (l,0). 

We will now summarize the above discussion in the following definition of 
the renormalization operator acting on generating functions originally due to the 
authors of (Eckmann et al 1984): 

Definition 2.1 

T^EKw[s](x,y) = p~ 1 s(z(x,y), Ay), where 

= s(Xx,z(x,y)) + s(Xy,z(x,y)), 

= a (A,l) +s(0,l) and p = 8^(1,0). (2.6) 

Definition 2.2 The Banach space of functions s(x,y) = J2Tj=o Ci j( x ~ T ) l (v ~ 
t) 3 , analytic on a bi-disk 

\x-t\< p, \y- t\ < p, 

for which the norm 

oo 
i,j=0 

is finite, will be referred to as A T (p). 

A T s (p) will denote its symmetric subs-pace {s G A T (p) : si(x,y) = s±(y,x)}. 
We will also use the shorthand notation 

A(p) = A°(p), M P ) = A° s (p)- 

As we have already mentioned, the following theorem has been proved with 
the help of a computer in (Eckmann et al 1982) and (Eckmann et al 1984), and 
later in (Gaidashev et al 2011). We will quote here the version of the theorem 
from (Gaidashev et al 2011), since we will use the same functional spaces as the 
proofs in that paper. 

Theorem 1 There exists a polynomial s a : C 2 i— > C, such that 

i) The operator TZekw is well-defined, analytic and compact in B t {sq) C 
A s {p), with 

p=1.75, r = l.lxl0~ 10 . 



ii) There exists a function s* G B r (so) C A s (p) such that 

TIekw[s*} = s*. 
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in) The linear operator D1Zekw[ s *} has two eigenvalues outside of the unit 
circle: 

8.72021484375 < 5 1 < 8.72216796875, 5 2 = t-, 

A* 

where 

-0.248875313689 < A* < -0.248886108398438. 



iv) The complement of these two eigenvalues in the spectrum is compactly con- 
tained in the unit disk. The largest eigenvalue in the unit disk is equal to 
A*, while 

spec(DTZ E Kw[s*])\{Si,5 2 ,X*} C {z G C : \z\ < 0.1258544921875}. 



In this paper we will adopt a normalization condition different from (j2.6|) 
however. Specifically, p will be defined from the normalization condition 

Kekw[s](0,0) = 1, 

i.c 

m[b] = *(*(0,0),0). (2.7) 

Definition 2.3 Define l T (p) and lj(p) to be the subsets of A T (p) and A T {p), 
respectively, of functions normalized in the following way 

l T (p) = {sEA T (p):s(0,0) = l}, 

Zl(p) = {seA T s (p):s(Q,0) = l}, 

J(p) = {s £ A(p) : s(0,0) = 1} , 

l s (p) = {s e A s {p) : s{0,0) = 1} . 

Our main result concerning the operator 12.11 will be the following 

Main Theorem 1 There exists a polynomial so € ^{p), p = 1.75, of degree 9 
such that 

i) the operator {jO]) is analytic in B s (s ) C A(p) with 5 = 0.00405550003051758; 

ii) for all s G B$(sq) with real Taylor coefficients, the scalings A = X[s] and 
p = p[s] satisfy 

0.000406771898269653 < p < 0.120654106140137, 
-0.276069164276123 < A < -0.222213745117188; 

in) the operator (jg.jp is compact, with7ZEKw[s\ £ Z(p')> p' = l-0699462890625p, 
for all s G Bs(sq). 
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Remark 2.4 We would like to note that the renormalization fixed point s* sat- 
isfies 

II s* - so Hp < .00368565320968628, 

(which is a rigorous bound), i.e., as expected, according to Theorem^ lies in the 
analyticity domain of renormalization. 

Remark 2.5 All the numbers quoted in the Theorem and the Remark above are 
representable on a computer. 

Although analyticity and compactness of the renormalization operator have 
been already proved in (Eckmann et al 1984) (together with a stronger result of 
existence of the fixed point), we would like to reiterate that the goal of the paper 
will be to obtain these results with a lighter machinery than that of (Eckmann 
et al 1984). We also obtain a better bound on the size of the neighborhood of 
analyticity of (|2.1|) . 

3 Interval arithmetics in R and IR 2 

We will now give a very brief summary of interval operations in K and M 2 . For a 
more complete treatise of "standard sets" and operations on them, an interested 
reader is referred to an excellent review (Koch et al 1996). 

A computer implementation of an arithmetic operation r\-j^r<i is +, — , * or 
/) on two real numbers does not generally yield an exact result. The "computer" 
result is a number representable in a standard IEEE floating point format (cf 
(IEEE 1985)). Such numbers are commonly referred to as "representable". In 
the 80-bit extended precision IEEE arithmetics, a number is represented with 80 
bits of memory: 1 bit for the sign of the number, 15 bits for an exponents, and 
64 bits for a mantissa. A real representable number r is of a the form 

r = ±l.m 2 e ~ (2l4 ~ 1} , 

where ± is chosen according to whether the sign bit is or 1, m is a base 2 
mantissa given by the sequence of 0's and l's associated with the 64 bits of the 
mantissa, and the exponent < e < 32767 is defined by the state of the 15 
exponent bits (the state of 15 ones is reserved to represent "overflows"). The 
representable number is given by the sequence of 80 zeros. The set of all 
representable numbers will be denoted 1Z. 

Now, let r\j^r<i be a mathematically legal, nonzero, arithmetic operation. 
The result of the true arithmetic operation r\j^r2 might not be a representable 
number. However, we can instruct the computer to attempt to round this oper- 
ation either to the nearest representable number, up or down (our choice was to 
always round up). This might not be possible: the result of rounding up might 
have the exponent e > 32767 (overflow), or e < (underflow). In both cases 
the computer is instructed to raise an exception, which is appropriately handled 
(either by terminating the program, or restarting with a different set of param- 
eters). If, however, the result of rounding up is representable, the output of the 
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computer implementation of the arithmetic operation is an upper bound on the 
true result of the operation. We will refer to such bound as "the upper bound" . 

We will define the standard sets in R to be the collection of all closed real 
intervals I[x,y] = {r G M.\x < r < y}: 

std(R) = {I[x, y] £R: x,y £K}. (3.8) 

If I[x,y], I[xi,yi] and I[x2,y2] are in std(R) then we can use the rounding 
up described above to obtain bounds on the arithmetic operation on these sets 
as follows: 

1) unary minus: —(I[x,y]) = I[—y, —x]; 

2) absolute value: \I[x,y]\ = I[l, r], where I = max{0, x, —y}, r = — min{0, x, —y}; 

3) addition: I[xi,y±] + I[x2,y2] = I[%3,y3], where 2/3 is the upper bound on 
Ui + 2/2) while — X3 is the upper bound on —x\ + (—X2); 

4) subtraction: I[xi,yi] - I[x 2 ,y2] = I[xi,m] + (-^[^2,2/2]); 

5) multiplication: I[xi,y±] • I[x2,y2] = ^[^3, 2/3] , where j/3 is the maximum of 
the upper bounds on x\ ■ X2, X\ ■ y2, 2/1 • x 2 and 2/1 • 2/2 ; while — X3 is the 
maximum of the upper bounds on (— x±) • X2, (— X\) ■ y2, (—2/1) " x 2 and 
(-2/1) • 2/2; 

6) inverse: if x\ ■ y\ > 0, then I[xi,yi] inverse = I[x2,y2], where 2/2 is the 
upper bound on 1/xi and — X2 is the upper bound on l/(— 2/1); 

7) division: if x 2 -2/2 > 0, then /[xi, j/i]/J[x 2 , 2/2] = ^[ari, 2/1] ■ (-^[^2, 2/2] inverse). 
A standard set in R 2 is, naturally, a direct product of two standard sets in R: 

std(R 2 ) = {/[xi,yi] x I[x 2 ,y2],x 1 ,X2,yi,y2 € 11}. (3.9) 

The arithmetic operation on these standard sets in R 2 are reducible to those 
on reals in an obvious way. 

To obtain a bound on a algebraic and transcendental function, one can use 
their Taylor series together with a bound on the remainder. In our proofs we will 
require only three such functions: exp, In and y/~. All of them are implemented 
via a finite truncation of their Taylor series with a bound on the remainder. 

4 Contraction Mapping Principle 

We will now outline a rather general method for finding a fixed point of a hyper- 
bolic operator in a Banach space via its approximate Newton map. 

Let C be an operator analytic and hyperbolic on some neighborhood M in 
a Banach space Z. Suppose, that one knows its approximate hyperbolic fixed 
point Z e M. Set 

M= [I-DClZo]}- 1 , 
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and for all z, such that Z$ + Mz € Af, 



N[z] =z + C[Z + Mz] - (Z + Mz). 

Notice, that if z* is a fixed point of N, then Zq + Mz* is a fixed point of C. 

The linear operator I — DC[Zq] is indeed invertible since DC is hyperbolic at 
^0- If Zq is a reasonably good approximation of the true fixed point of C, then 
the operator N is expected to be a strong contraction in a neighborhood of 0: 

DN[z] = I + DC[Z + Mz}- M - M 

= [M- 1 + DC[Z + Mz] — I] ■ M 

= [I - L>C[Z ] + £>C[Z + Mz] — I] ■ M 

= [DC[Z + Mz] - £>C[Z ]] • M. 

The last expression is typically small in a small neighborhood of 0, if the norm 
of M is not too large (if M is large, one might have to find a better approximation 
Zq and take a smaller neighborhood of 0). The following well-known Theorem 
specifies a sufficient condition for existence of the fixed point: 

Theorem 4.1 (Contraction Mapping Principle) 

Suppose that the operator N is well-defined and analytic as a map fromN C Z 
to Z, where Z is some Banach space. Let Zq £ M and B$(Zo) C N (an open 
ball of radius 8 around Zq ) be such that 

\\DN[Z]\\ <V<1, 

for any Z G Bs(Zo), and 

\\N[Zo]-Z \\<€. 

If e < (1 — V)S then the operator N has a fixed point Z* in Bs(Zq), such that 



We will use the above Contraction Mapping Principle in several instances in 
our proofs below. In all those cases we will be verifying the hypothesis of the 
Contraction Mapping Principle using mild computer assistance. 



5 Analyticity of renormalization 

We consider the space A(p), p = 1.75. Define 

8 (x, y) = D (x)y 3 + A (x)y 2 + B (x)y + C (x), 
a polynomial of degree 3, with 

6 6 6 6 

D (x) = dix\ A (x) = ^2 ai x\ B (x) = hx\ C (x) = ax\ 

i=0 i=0 1=0 i=0 

where the coefficients of these polynomials are as follows (the numbers given 
in the table are highly accurate approximations of the representable numbers 
actually used in the programs): 
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i 


Ci 







1.00000000000000000 


-2.42962369607899157 x 10" 1 


1 


-1.02761956458970711 


5.87327440047455615 x Wr' 1 


2 


2.93663720023727808 x 10" 2 


-5.93710236103475834 x 10" 3 


3 


rj — 

-1.87658664952086400 x 10~ J 


6.09332694202817819 x 10~ 4 


4 


1.40668294317213841 x 10" 4 


p 

-6.46957663100331420 x 10~ 5 


5 


F 

-1.18664608613747513 x 10" 5 


T> 

7.02844653606969302 x 10~ fa 


6 


1.06935654680404746 x 10~ e 


-7.75814237637266867 x 10~ 7 


i 


di 


di 





-8.77647505670140721 x 10" 1 


-1.46791670728014469 x 10~ 3 


1 


-5.62975994856259201 x 10~ 3 


5.62673177268855366 x 10" 4 


2 


9.13999041304226728 x 10~ 4 


-1.29391532620066284 x 10" 4 


3 


-1.33251463600020727 x 10~ 4 


2.45371180652649971 x 10" 5 


4 


1.84028385489487478 x 10" 5 


-4.18641657324405651 x 10~ e 


5 


-2.46417657757601024 x 10~ B 


6.68319468409332288 x 10" 7 


6 


3.23459112039180563 x 10~ 7 


-1.01939399249366523 x 10~ 7 . 



A numerical study demonstrates that so is a relatively good approximation 
of the renormalization fixed point in A s (p) for a rather wide range of radii p 
(p « 0.8 . . . 2.5, below p will be fixed to be 1.75). 

In this Section we will prove Theorem [TJ This Theorem is proved with a 
mild aid of the computer. The computer assistance is restricted only to interval 
arithmetics on real numbers, (which, given sufficient time, in principle, can be 
done "by hand"). 

We will start by noticing that a solution of the equation 

s(x,Z') + S (y,Z')=0 

with s = so is given by 



(A (x) + A (y)(Co(x) + C Q (y) 
(B (x) + B (y)y 

where Cq{x) is a bound on Cq{x) + y 3 D(x), specifically: 

C (x) = C (x) + CD Q (x), C e B p3 (0). (5.10) 

In the formula above we choose a branch of such that the range of Z Q is 
expected to be contained in D p . 

We would like to remark that we have intentionally chosen a polynomial so of 
a sufficiently high degree so that it would be close to the fixed point s* found in 
(Gaidashev et al 2011). This will later imply that the fixed point s* is indeed in 
the domain of analyticity and compactness of IZekw- At the same time, taking 
a polynomial of degree higher than 3 in y is superfluous for the purposes of the 



B (x) + B (y) 
Z °^ y) ~ 2(A (x) + Mv)) 
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computation of Z' . Therefore, these extra degrees (polynomial Dq) have been 
included as an "error term" in (|5,10p . It is clear, that 

C (x) + y 3 D (x) GCoW 

for all (x,y) such that \x\ < p, \y\ < p. 

We shall now check that for the given values of coefficients the function Z' is 
analytic on a bi-disk of radius r = 0.483119964599609 centered at point (0, 0). Be- 
low Z(r) will denote the Banach space of functions analytic on {(x, y) : \x\ < r, \y\ < r} 
equipped with the sup- norm, denoted by | • | r . We have the following 

Lemma 5.1 Zq is in Z(r) for r = 0.483119964599609, and satisfies 

\Z'v\ r < 1.562789916992188. 

Proof. First, define constants vo, (3q and 70, and quartic polynomials v, (3 
and 7, such that v(0, 0) = /3(0, 0) = 7(0, 0) = by setting 



v + v(x,y) 
70 + 7(x,y) 



A (x) + A (y), 
B (x)+B (y), 
C (x) + C (y). 



Then 

z 'o(x,y) 
where 



Po + P(x,y) 
2(v + v(x,y)) 



i^- + F 1 (x,y) - F 2 (x,y) - 1 



(5.11) 



Fi(x,y) 
F 2 (X,Y) 



4v ~fo 2(3 (3(x,y) + f3 2 (x,y) 
Pi ^3l + 2f3of3(x,y)+^3Z(x,y) , 
v j(x, y) + v(x, y)7o + v(x, y)j{x, y) 



/3 2 + 2/3 /3(x,y)+/3 2 (x,y) 
Norms of both F\ and F 2 are elementary: 



4|^7o| |2/3 /3 + /3 2 | r 
Pi /3 2 -|2/3 /? + /3 2 



1 2|r " /? 2 -|2/3 /3 + /3| 2 



Notice that all functions here whose norm has to be evaluated are explicit poly- 
nomials. We have estimated the above norms using the interval arithmetics on a 
computer. Next, we use these bounds in the following expression 



\Z \ r < 



IA)| + 



— j— - max 

V\r) i=±l 



2(K 

and evaluate it on a computer 



4— ^- + 3\Fi\ r +j\F 2 \ 



1 



1 
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□ 

Below, we will use the following shorthand notation for the supremum of Z' Q : 

t = \Z' \ r . (5.12) 

Notice that if one writes 

s(x,y) = a{x,y) - ^y-, 

where r is chosen so that it has an inverse branch v with a range in D p (it 
is sufficient to choose r a quadratic function), then a solution of equation (0) 
satisfies 

Z'{x, y) = u[a(x, Z'(x, y)) + a(y, Z'(x, y))]. 

(this "trick" has been shown to me by Hans Koch). 

We will now specify the choices of r and a for sq- Set 

s (x,y) = a {x,y) 

<?o(x,y) = a(x)y 2 + b(x)y + c(x), and S (x, y) = a (x + p, y), 

T (.y) = -2a y 2 - 2b y, 

-b + \/b\- 2a x 

v(x) = . 

2a 

It can be readily verified that Z Q (x,y) = z^[Ilo(x, Z (x, y)) + Y>o(i), Z'^x^y))]. 
We will next demonstrate that the operator Ch, defined on Z(r) by setting 

c h [z'\ = v [s o (n 1; z') + ho (p u z') + s o (n 2 , z') + ho (p 2 , z')] , 

where 

Pi(x,y) = x +p, and P 2 (x,y)=y+p, 
has a fixed point for sufficiently small /i's. 

Proposition 5.2 Let 5, p and r be as in Theorem [I] and Lemma \ 5.1\ Then, 
for all s G Bs(sq) C A(p) the operator C S _ S( , has a unique fixed point Z' s in 
B e (Z ) C Z{r) with e = 0.01465, and the map s h-» Z' s is analytic from Bs(so) to 
B e (Z' ). 

Proof. Define an operator 

N h [z] =z + C h [Z + Mz) - (Z' + Mz), M = [l - DC [Z }} ^ . 

We will demonstrate that this "Newton map" has a fixed point in a neighborhood 
of 0. 
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We will first estimate the norm of N h [0] = C h [Z' ] - Z : 



\Ch[Z' Q ] - Z' \ 



\v [s o (ni, z' Q ) + h o (iii, z£) + s o (n 2 , z' Q ) + ho (n 2 , z' Q )] 
- ^[s (ni,z ) + £o(n 2 ,z )]| r 

2^ - 2 «o [s (ni,^) + ^(Pi,^) + s (n 2 ,z ) + fr(P 2 ,z )] 



- y/%-2ao [SoOl^ + SoCn,,^)] 



2a 



|P(l)-F(0)| r , 



where 

TO 



&3 



2a [E (ni, Z ) + (h(P 1 , Z' Q ) + S (n 2 , Z ) + (h(P 2 , Z' )\ 



is an analytic function from I$r to Z(r) (here and below, B# stands for an open 
disk of radius 22 in C), with 



22 < 



b 2 Q - 4a # - 2|a ||e (IIi,Z£)| r - 2|a ||G (n 2 , Z' ) 



(5.13) 



2|ao||/i(Pi,^)| r + 2\a \\h(P 2 ,Z^)\ r 
where 

9 = S (0, 0), 9 (x, y) = S (x, y) - 9 . 

Below, we will denote the (positive) real number 6q — 4ao#o by c. A straightfor- 
ward Cauchy estimate yields: 

4\a \5\/2c 



\F(1)-F(0)\< 



22-1 



We will denote 



S |i?(C)| - c-4|ao||e |(r )t )-4|ao|<5" 
2SV2c 



e = 



(5.14) 



c — 4|ao||©o|(r,t) — 4|ao|<5" 
This expression has been shown to satisfy 

e < 0.0137615203857422. 

As a next step we will provide a bound on the derivative of the operator N^: 

\DN h [z}\ = \[DC h [Z' + Mz]-DC [Z , ]]M\ 
< \[DC h [Z' + Mz]-DC [Z' }]\\M\. 

To estimate the norm of M we first verify that |DCo[Z ]| < 1: 

d 2 & (u 1 ,z^ + d 2 e (n 2 ,z^ 



DC, 



< 



2^c - 2a [6 (ni, Z' ) + 9 (n 2 , Z' Q )\ 
\d2@o\ {r , t) 



yjc-4\a \\e \ ( 



r,t) 
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A computer-aided evaluation of the last expression indeed shows that it is less 
than 1. This implies that M can be found as the limit of a convergent series 
whose norm is easily bounded: 



\M\<- 

Next, 

DC h [Z']-DC [Z' ] = 



< 1.0430755615234375 = M. 



(5.15) 



+ 



\DC [Z' ]\ 

^e (ni, z 1 ) + a 2 e (n 2 , z') + d 2 h(p u z') + d 2 h{p 2 , z') 
2y/ c - 2a [e (ni, z>) + e (n 2 , z>) + h(p u z>) + h(p 2 , z<)\ 

^eo(ni,z^) + ^e (n 2 ,z^) 
2^ c - 2a [e (ni, zl) + e (n 2 , z )] 

^e (ni, zp + a 2 e (n 2 , z') + d 2 /i(Pi, go + d 2 h(p 2 , z') 
2^/c - 2a [e (iii, z>) + e (n 2 , z 7 ) + h(p 1 , z>) + /i(p 2 , z>)\ 

d 2 Q Q {Ii 1: Z') + d 2 Q Q {Ii 2 ,Z') 

2y/ c - 2a [e (ni, z 1 ) + e (n 2 , z 7 ) + h{p u z>) + /i(p 2 , z')] 
,9260(11!, z') + g 2 e (n 2 ,z') 

2^/c - 2a [6 (ni, Z') + 9 (n 2 , Z') + h(P 1} Z>) + /i(P 2 , Z')] 

d 2 e (u 1 ,z') + d 2 @ (u 2 ,z') 



+ 



2y/ c - 2a [e (ni, z') + e (n 2 , z')] 

freoCn^zQ + 3260(112,^) 
2^c - 2a [e (ni, z') + e (n 2 , z')] 

^eo(ni,z^) + ^eo(n 2 ,z^) 
2y/ c - 2a [e (ni, z ) + e (n 2 , z )]_ 
= /i + j 2 + j 3 . 

We will estimate norms of the three expressions in brackets separately. 

|<9 2 /i|( rjS ) mS 



(5.16) 



|/l|r < 



where 



^c-4|ao||e |( r ,a) -4|a ||/i|( r , fl ) y /c ~ 4 l a oll o|(r,s) -4|a |<5 



(5.17) 



s = t + |M|e, m = n- 



n— 1 



and n 



MS] 



We use a Cauchy estimate for the function 



1 



- 2a [e (ni, zo + e (n 2 , z') + CKPi,z') + c/*(p 2 , ^0] 

to bound |/ 2 | r . Notice, that, F 2 is also analytic on Or, but gets unbounded as £ 
approaches the boundary of this disk. Therefore, to bound \I 2 \ r , we use a Cauchy 
bound for F 2 on a smaller disk of radius 
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c- 4|a 1 |©o|( r , s ) 1 c - 4|a | |6 |( riS) - 4|a |(5 
4\a \5 3 4|a |<5 



1 4|ao|<5 |<92©oL ,i 

\h\r<\d2&0\ (r , s) -~— SUp |F 2 (C)|< " 



R ~ 1 \^ R 2 (| (c - 4|a Q | |O | (r , s) - 4|a |^ V 

Finally, to estimate |i3| r , we use a Cauchy bound for the function 

%e (iii, z' + (Mz) + <9 2 e (n 2 , z' + (Mz) 



F 3 (C) 



2y/c - 2a [6 (ni, Z' Q + (Mz) + G (n 2 , Z' Q + (Mz)} ' 
analytic on with 

\M\e ■ 

\h\ r = \F 3 (1) - F 3 (0)| < J— sup |F 3 (C)| < \Mh^^k^=. (5.18) 

^ _1 ICI<A p s yc-4|a ||©o|(r lP ) 

Individual norms |M|, |/i| r , |/ 2 | r and |/3| r and their sum have been estimated 
on a computer to produce: 

\DN h (z)\ < 0.0125999450683594 = V (5.19) 

for all \\h\\ p < 5 and \z\ r < e. 

Estimates (|5.14p and ()5. 19|) demonstrate that 

■ <(!-!>)£. 

and the claim follows from the Contraction Mapping Principle. 
□ 

At the next step we obtain bounds on the scaling parameters A and [i. 
Recall that the scaling A solves equation (|2.6p . 
The prove Theorem Q] we will require the following 

Lemma 5.3 Equation \2.b\ has a solution A = A [a] for all s G B§(so), and 
satisfies 

-0.276069164276123 < A < -0.222213745117188 

for all s G Bs(sq) with real Taylor coefficients. 

Furthermore, the map s t— >■ X[s] is analytic on Bs(so). 

Proof. Write, 

si (x,y) = Ai (y)x 2 + B t (y)x + C x (y) , 
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where 

M (y) = c 2 + b 2 y + a 2 y 2 + d 2 y 3 , 
B\ (y) = ci + biy + a 1 y 2 + d x y 3 , 
C\ (y) = c + b y + a y 2 + d y 3 . 

Denote 

Notice, that the equation ()2.6p with s = s\ can be solved exactly (a quadratic 
equation for A). Furthermore, suppose s = s\ + h, \\h\\ p < \\so — si\\ p + 5 = S, 
then X[s] solves 

si(A[s], 1) + si(0, 1) = -h(X[s], 1) - h(0, 1). (5.20) 

Write 

A[s] = A[si] + AA, 

then AA is a solution of 

A^l) ■ (AA) 2 + [2Ai(l) • A[si] + 5i(l)] • AA + C = (5.21) 

with C € © 2( 5 (an open disk of radius 25 in C). For all s with real coefficients, 
(|5.2ip is a quadratic equation for AA with real coefficients, whose solution gives 
the required bound. 

□ 

The last ingredient in the proof of analyticity of renormalization is a bound 
on the scaling [i. 

Lemma 5.4 Let /x[s] be as in equation (2.1). Then the map s i— /i[s] is analytic 
on B$(so), and /i[s] satisfies 

0.000406771898269653 < n < 0.120654106140137, 

for s G Bs(so) with real Taylor coefficients. 

Proof. The bound is straightforward: 

fi[s] = 8(^(0,0), 0) 

€ s (^(0,0),0)+M^(0,0),0). 

Now, it is clear, that whenever s has real coefficients, 

li[s] G s (Z' (0, 0) + I[-\M\e, \M\e],0) + I [-5, 5} , 

where, as before, I[x, y] stands for a closed real interval with endpoints x and y, 
while the rest of the quantities are as in Lemma 15.31 and Prop. 15.21 

□ 

Now, part i) of Theorem [T] follows immediately from Proposition 15.21 and 
Lemmas 15.31 and 15. 
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6 Compactness of renormalization 

We will now outline the proof of compactness of the renormalization operator. 

Proposition 6.1 For every s £ B$(so) C A(p) its renormalization TZekw[$] is 
in X{np) with k = 1.0699996948242188.. 

Proof. Let r' = nr and p' = up. Verification that |A[s]|p'| < p is straightfor- 
ward. The "difficult" part is to demonstrate that \Z' s \ r / < p. To show this, we 
will use the fact that Z' s = lim^oo iV s _ so [0], 

To that end we first verify that Z' £ Z{r'). This amounts to checking on a 
computer that 

|uo| > \v\ r > and 1 -4^f > |Fi| r / + |F 2 | r / (cf. (fTTTTH ). 

Po 

Below we will use the following shorthand notation: t' = \Z'q\ t i . 
Set zq = 0, and define for all n > 

z n+1 = N h [z n ], Az n = z n+1 -z n , Z' n = Z' Q + Mz n 

(we have suppressed the dependence on s in Z' n for notational convenience) , then 

|Az | r / = \Ch[Z ] — Z \ r > 



and 



25v2c , „ . _. 

< ,i Mn | ( cf - equation (EH} , 



\Az n \ r > = \N h [z n ] - Nh[Zn-i]\ r ' 

= |Az n _! + C h \Z' n _ x + MAz n _i] - C h [Z' n _ x } - MAz n - X \ r , 

< | [C h [Z' n _ x + MAz n -i] - C h [Z' n _ x ] - DC h [Z' n _ x ]MAz n ^ x ] \ r , + 

| [DC h [Z' n _ x ] - DC h [Z^}} MAz n ^\ r , 
= T\ + T 2 . 

To evaluate T 2 we will require estimates similar to (|5.15p - (|5.18p . In fact, esti- 
mates (|5.15p ~ (|5.18p can be recycled after one substitutes r' for r, e' for e and 
s' = t' + \M\e' for s (here e' is an a priori bound on |.z n | r ', n > 0, again, veri- 
fiable a posteriori). We will not repeat the details, but rather use the following 
symbolic notation 

|T 2 | r / < C 2 |Az n _i| r ', 

keeping in mind that C 2 is computable on a machine. 

To estimate T x , we use a Cauchy bound for the following function: 



F(C) = ^c-2a [Q (n 1 ,Z' n _ 1 +(MAz n ^ 1 ) + h(P 1 ,Z' n „ 1 +(MAz n _ 1 ) + 
6 (n 2 , Z' n _ 1 + QMAz n ^ l ) + h(P 2 , Z^+CMAzn^)] , 
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analytic on Br with 

p-K-il 



R < 



\M\\Az n ^\ r r 
This results in 

" C V-^-|Ml|ALi|r'' (6 - 22) 



Finally, we have 



\Az n i2 



\Az n \ r > < Ci ' j— + C 2 |Az„_i| r / 

p-s' - \M\\Az n - 1 \ r > 

= K n -i\Az n -i\ r >\, 

K - r l Az »-l|r' ~ 
= Oi j h0 2 

p-s' - \M\\Az n -i\ r > 

(we check that < 1). 

Clearly, K„ is a monotone increasing function of |Az n | r ', therefore K n < 
K n -\. We can now verify that the sum 



^ = Y^\Az k [, 



k=0 

is convergent: 



S < ^if fc _i|Az fe _i| r / + |Az |r' 

fc=i 

oo 

< ^ Kfc_iKfc_ 2 | A^_ 2 | r ' + K \Az \ r < + |Az | 



fc=2 

oo 

< ^K fe |A. |, 

fc=0 

< [AgpLy 

Therefore, for all n > 



1 - A 



Finally, we verify that 



t' + y < p. 



□ 
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